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ABSTRACT 

The ionising stellar populations of eleven Hu regions in the spiral galaxies: 
NGC 628, NGC 925, NGC 1232 and NGC 1637, all of them reported to have 
solar or oversolar abundances according to empirical calibrations, have been 
analysed using stellar population synthesis models. 

Four of the observed regions in the sample, show features which indicate 
the presence of a population of Wolf-Rayet (WR) stars with ages between 2.3 
and 4.1 Myr. This population is sufficient to explain the emission line spec- 
trum of the low metallicity region H13 in NGC 628, taking into account the 
uncertainties involved in both observations and model computations. This is 
not the case for the rest of the regions for which a second ionising popula- 
tion is required to simultaneously reproduce both the WR features and the 
emission line spectrum. Composite populations are also found for half of the 
regions without WR features, although in this case, the result is based only 
on the emission line spectrum analysis. 

For two of the regions showing WR features, no consistent solution is 
found, since the population containing WR stars produces a spectral energy 
distribution which is too hard to explain the emission of the gas. Several 
solutions are proposed to solve this problem. 

Key words: galaxies: abundances - galaxies: HII regions - galaxies: stellar 
content - stars: WR stars 
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1 INTRODUCTION 

Metallicity affects stellar evolution in at least two different ways: by increasing the opacity of 
the stellar material and through the strengthening of the wind driven mass loss in high mass 
stars. The former implies a lower effective temperature for the atmospheres of ionising stars 
of higher metal content even at zero age, while the latter can severely affect the evolution 
of the most massive stars leading, in the most extreme cases, to the almost complete loss of 
their outer envelopes. 

The first of these effects should be readily observable as an inverse correlation between 
metallicity and temperature of the ionising radiation for HII regions of different metallicities. 
In fact, there is a general agreement about the hardening of the ionising radiation in regions 
of low metal content (e.g. Campbell et al. 1986; Cervino and Mas-Hesse 1994). The inverse 
situation, however, i.e. the softening of the ionising radiation in regions of high metal content, 
is more difficult to establish since, in these cases, the determination of both metal abundance 
and ionising temperature is rather difficult. For the first one we have to rely , in general, 
on uncertain semi-empirical calibrations when the involved abundances are higher than 
about 1/2 solar. For the second one, there are no suitable nebular emission line diagnostics. 
Recently, Bressolin, Kennicut & Garnett (1998) have investigated a possible relation between 
ionising stellar temperature and metallicity in giant HII regions in spiral galaxies through 
the use of the softness parameter r\ (Vilchez & Pagel 1988) concluding that their results are 
consistent with a significant decrease in mean stellar temperatures of the ionising stars with 
increasing metallicity. This decrease is however difficult to quantify since in this work the 
metallicity is derived from the abundance parameter R 2 3 (Pagel et al. 1979) and both 77' and 
R23 depend to some extent on the degree of ionisation of the nebula (Diaz et al. 1991). On 
the other hand, the method can be applied to a large number of HII regions and therefore 
can provide trends statistically significant. 

A different approach consists in deriving the metal abundance from the direct determina- 
tion of the electron temperature in the nebula and the ionising temperature from the fitting 
of single star photoionisation models. The method requires both high quality spectroscopy 
and detailed modelling. In a previous paper (Diaz et al. 2000a; DCTG00) we have applied 
this procedure to HII regions in the spiral galaxy NGC 4258. The combination of these 
data with similar ones in the literature also seems to point to a decrease of mean ionising 
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temperature with increasing metallicity. Surprisingly, this is so even in the presence of WR 
stars. 

The presence of these stars provides an additional constraint to characterise the ionising 
stellar population. In the last years, stellar evolution models for massive stars taking into 
account mass loss at different metallicities have been calculated in order to reproduce the 
observed features of WR stars. The model predictions are however quite different depending 
on the assumed mass loss. In fact, models that assume a standard mass loss rate underpredict 
the number of galactic WR stars in comparison with observations (Maeder & Meynet 1994) 
while an enhanced mass loss rate fits most WR population properties except for the mass 
loss rate itself (Leitherer, Chapman & Koribalski 1997). Therefore, the detailed observation 
and modelling of giant HII regions showing WR features can help to shed some light on this 
very important issue. 

It would also be desirable to check the hardening of the ionising radiation due to the 
presence of Wolf-Rayet stars predicted on theoretical grounds (Perez 1997). These stars are 
supposed to change drastically the spectral energy distribution of an ionising star cluster at 
energies higher than ~ 40 eV. Models predict that this change depends both on the age and 
the metallicity of the burst (Schaerer & Vacca 1998; SV98) and, as a result of this, ionic 
elemental ratios are expected to change as well as ion-weighted temperatures. 

A puzzling result in DCTG00 was the relatively low effective temperatures found in 
the two ionising clusters with WR features, suggesting either that the rate of high energy 
photons is lower than expected or that the opacity in WR envelopes is highly efficient. This 
matter can be investigated through the observation of WR features in high metallicity HII 
regions where the effect of the mass loss rate on the effective temperature of the ionising 
stars is supposed to be most enhanced. 

In order to analyse the effects of metallicity on both stellar evolution and ionising ra- 
diation, we have analysed a sample of reported high metallicity Hn regions by combining 
evolutionary synthesis models for ionising populations and photoionisation models in a self- 
consistent way. In this paper a detailed analysis of the ionising populations of the observed 
regions is presented. 
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2 PHYSICAL PROPERTIES OF THE OBSERVED HII REGIONS 

The observations of the sample HII regions analysed here: 4 in NGC 628, 2 in NGC 925, 4 
in NGC 1232 and 1 in NGC 1637 are described in detail in Castellanos, Diaz & Terlevich 
(2002; Paper I) and consist of spectrophotometry in the 4000-9700 A range. These data 
were complemented with those by Van Zee et al. (1998) in order to provide a complete 
set of emission line ratios in the spectral range from 3500 to 9700 A. For every region the 
physical conditions of the ionised gas- particle density, electron temperatures, ionic and 
total abundances - and functional parameters - ionisation parameter and temperature of 
the ionising radiation - as derived from the application of single star photoionisation models 
were determined and are summarized in Table 1 for the regions analysed here. 

For each of them we have calculated the Ha luminosity, L(Ha), the number of hydrogen 
ionising photons, Q(H), the angular effective diameter, 0, the filling factor, e, the mass of 
ionised hydrogen, M(HII), and the mass content in the ionising clusters, M*. These quan- 
tities have been derived according to the expressions given in Diaz (1998), for our adopted 
distances: 7.3 Mpc for NGC 628, 8.6 Mpc for NGC 925, 21.5 Mpc for NGC 1232 and 7.8 
Mpc for NGC 1637 (see Paper I) and are listed in Table 2. 

The Ha fluxes for regions H13, H3, and H5, uncorrected for reddening, are 669, 95 and 
62 x 10~ 16 erg s _1 cm -2 respectively, a factor of about three lower than those measured by 
Kennicutt & Hodge (1980) from Ha photometry (2089, 372 and 145 x 10 -16 erg s" 1 cuT 2 ). 
This factor reduces to about 2.4 if we take into account that their Ha filter also includes 
the [Nil] lines which amount to about 1/3 of Ha. 

Ha photographs (Hodge 1976) reveal these regions to be rather extended, and considering 
that our observations have been obtained through a 1'.'03 width slit, this comparison implies 
that our derived values should be considered as lower limits. At any rate, our value for H13 
agrees well with that measured by Van Zee et al. (Van Zee, private communication) which 
was obtained through a 2 '.'0 diameter aperture. 

For the rest of the regions we have not found any published photometric data. 

The angular diameter of the emitting regions as well as the filling factors can be obtained 
from the ionisation parameter, u, the reddening corrected Ha flux, F(Ha) and the derived 
electron density, n#, under the assumption of spherical geometry. Then 

u = Q(H)/47icR 2 n H 
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where Q(H) is the number of hydrogen ionising photons per second, R is the radius of the 
ionised region and c is the speed of light. 



The corresponding angular diameter is then found to be: 
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where F(Ha) is in units of erg cm _2 s _1 , n# is in units of cm -3 and <fi is in arcsec. For 

the regions for which only upper limits to the electron density could be derived, a value of 

n#=40 cm" 3 has been assumed. Angular diameters thus calculated, do not depend on the 

assumed distance to the galaxy. For all the regions, we obtain angular diameters between 

0.6 and 2 arcsec. 

Filling factors for each observed region can also be determined from the reddening cor- 
rected Ha flux and the derived ionisation parameter, according to the expresion: 

' F{Ha) \ /UhV 

v io- 14 J uov 

where D is in Mpc and aB{H°,T) is the case B recombination coefficient for hydrogen for 
which we have used a value of 3.76 x 10~ 13 cm 3 s _1 , corresponding to T= 7000 K and n#= 
100 cm -3 (Osterbrock 1989). Filling factors thus derived depend inversely on the galaxy 
distance, D. Values between 0.008 and 0.13 are found for our adopted distances. These 
values are very similar to those derived for the observed HII regions in NGC 4258 (see 
DCTG00) and M 51 (Diaz et al. 1991). 

The corresponding masses of ionised hydrogen range from 1.5 x 10 4 M for region H4 
in NGC 628 to 8.8 x 10 4 M for region CDT4 in NGC 1232. 

Regions CDT1, CDT3 and CDT4 in NGC 1232 have Ha luminosities greater than 10 39 
erg s -1 and hence can be classified as supergiant HII regions as defined by Kennicutt (1983). 
The rest of the regions have Ha luminosities typical of HII regions in early spiral galaxies, 
although all of them are greater than 10 37 erg s _1 requiring more than a single star for their 
ionisation (Panagia 1973). 

For our observed regions, all values for Q(H) are greater than 10 49 photon s -1 . For regions 
CDT1, CDT3 and CDT4 in NGC 1232 values as high as 10 51 photon s" 1 are found. In the 
absence of dust and if no photons are scaping the region, a lower limit for the mass of 
the ionising clusters can be estimated by means of the H/3 measured equivalent width and 
the Ha luminosity for each region (see Diaz 1998). Assuming a Salpeter IMF with upper 
and lower mass limits of 100 and 0.8 M respectively, the estimated ionising cluster masses 
range from 1.7 x 10 3 M for region H4 in NGC 628 to 1.80 x 10 5 M for region CDT1 in 

© 2000 RAS, MNRAS 000, §-?? 



6 M. Castellanos et at 

NGC 1232. It can be inferred from these results that all the regions are ionised by relatively 
small clusters, except for the three supergiant HII regions in NGC 1232. Regions H3, H4 
and H5 show a ten to one ratio between their ionised hydrogen mass and the mass of the 
ionising stellar clusters which may be then very young. Their high observed EW(H/3) values 
of about 200 A support this claim. 

If a given fraction of ionising photons were absorbed by dust grains before they reach 
the neutral hydrogen, then the number of Lyman continuum photons emitted Q(H) em by 
the star cluster would be higher : Q(H) = f .Q(H) em with f < 1. Therefore, our assumption 
of f = 1, i. e absence of dust, implies that Q(H) em and correspondingly M* are lower limits. 
The rest of the derived nebular parameters: angular diameter, filling factor and mass of 
ionised hydrogen depend on the extinction corrected Ha flux, the ionization parameter and 
the hydrogen density of the emitting gas. All these quantities being observables, they are 
not affected by the presence of absorbing dust. If f is less than unity, Q(H) em would be larger 
than the derived Q(H) and obviously if all these emitted photons were ionising the gas we 
should observe a much higher Ha flux and the nebular parameters would be different. 

2.1 Wolf-Rayet populations 

Relatively prominent Wolf-Rayet features have been identified in several of the observed HII 
regions: H13 in NGC 628, CDT3 and, at a fainter level, CDT1 and CDT4 in NGC 1232. In 
all of them, the typical blue 'bump' around A 4660 A is observed and a red 'bump' around A 
5800 A is also present in regions H13 (NGC 628) and CDT3 (NGC 1232). Details concerning 
these WR features can be found in Paper I and are summarized in Table 3. 

2.1.1 Region HIS in NGC 628 

Assuming a distance to NGC 628 of 7.3 Mpc (Sharina et al. 1996), and a constant extinc- 
tion value through this region of c(H/3) = 0.29, as derived from the Balmer and Paschen 
recombination lines, the total luminosities of the blue and red bumps are (1.8 ± 0.2) x 10 37 
erg s -1 and (6.5 ± 0.6) x 10 36 erg s -1 respectively. The former value comprises the features 
of Nm AA 4634, 4640, and Hen A 4686 A lines. The contribution of the Nin lines to the blue 
WR bump is metallicity-dependent according to Smith (1991). In our case this contribution 
represents 0.45 times the total emission. Hence, the derived value for the stellar Hell line 
luminosity is: 
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L(HeII A4686) = (9.5 ± 0.4) x 10 36 ergs' 1 

Based on the characteristics of the WR features we can conclude that the observed WR 
stars can be classified as WN7 and using the calibration of Vacca & Conti (1992), 6 WN7 
stars are found in region H13. 

2.1.2 Region CDT3 in NGC 1232 

In this case, the total luminosities of the blue and red bumps, assuming the distance to 
NGC 1232 of 21.5 Mpc adopted by Van Zee et al. (1998) and our derived reddening constant 
c(H/3) = 0.32, are (4.3 ± 1.0) x 10 37 erg s" 1 and (9.4 ± 2.5) x 10 36 erg s" 1 respectively. 
The derived value for the stellar Hell line luminosity is: 

L(H el 1X4686) = (1.95 ± 0.50) x 10 37 'ergs' 1 

We classified the dominant WR stars as of WN8 type. Therefore using the calibration of 
Vacca & Conti (1992), 12 WN8 stars are found in this region. The detection of CII, CIII at 
A 5140A could be explained by the presence of early or intermediate WC stars. It is only 
possible to give a lower limit to the number of WC stars because it is difficult to establish the 
dominant WC type (SV98). Hence, we have adopted a mean value for the average observed 
luminosity of WC4 stars at A 5808A of 3.0 x 10 36 erg s" 1 . Therefore, at least 3 WC stars 
could be found in this region. 

2.1.3 Rest of the regions in NGC 1232 

Fainter Wolf-Rayet features are found in regions CDT1 and CDT4 in NGC 1232. Region 
CDT1 shows a stellar Hell line luminosity of (9.5 ± 1.0) x 10 36 erg s _1 compatible with the 
presence of 6 WN8 stars. 

For region CDT4, the measured Hell line luminosity is (5.0 ± 0.8) x 10 37 erg s _1 . 
Accordingly, 30 WN8 stars would be necessary to explain the observed luminosity. 

3 MODELLING THE IONISING POPULATIONS 

The time evolution of the ionising energy output of a young star cluster can be followed 
with the help of stellar evolutionary tracks combined with adequate atmosphere models and 
they must incorporate the most massive stars. This has been made in the last years by 
several independent groups (e.g. Cervino & Mas-Hesse 1994; Garcia- Vargas, Bressan & Diaz 

© 2000 RAS, MNRAS 000, 0-?? 



8 M. Castellanos et at 

1995; Leitherer & Heckman 1995). In these so called "evolutionary synthesis models" the 
ionising stellar population is represented by a cluster of coeval stars formed with a given 
IMF which evolve in time along to theoretical evolutionary tracks in the H-R diagram and 
that is able to ionise the surrounding gas. Then, the combination of the individual stellar 
spectral energy distributions (SED) along a given isochrone representing a given cluster age 
is computed. These individual SED can be taken directly from observations or from stellar 
atmosphere models. Although the method is rather complex and the assumption of different 
stellar evolution and atmosphere models can lead to different predictions for the SED of a 
given cluster, there are several important conclusions which are common to the work by all 
groups mentioned above. In particular, a hardening of the ionising spectrum is predicted to 
occur between 3 and 4 Myr after the formation of the cluster as a result of the appearance 
of WR stars. In fact, the presence of these stars, which can be detected by the conspicuous 
high ionisation emission lines produced in their powerful winds (see section 2 above) provides 
important clues to constrain the age of the ionising population. 

SV98 have presented very detailed models of the WR population in young star clusters, 
at different metallicities from Z=0.001 (1/20 Z ) to Z=0.04 (2 Z ). Their clusters are formed 
according to a Salperter IMF with upper and lower mass limits of 120 M and 0.8 M re- 
spectively. They use the stellar evolution models by Meynet et al. (1994) which assume a 
mild overshooting and enhanced mass loss rate. These models have been shown to reproduce 
the observed WR/O star ratios in a variety of regions (Maeder & Meynet 1994). Regard- 
ing the energy output of main sequence stars they use the CoStar models (Schaerer & de 
Koter 1997) which include non-LTE effects, line blanketing and stellar winds, for stars with 
initial masses larger than 20 M and Kurucz (1992) plane-paralel LTE models, including 
line blanketing effects, for less massive stars. The atmospheres of evolved stars in the WR 
phase correspond to the spherically expanding, non-LTE, unblanketed models by Schmutz, 
Leitherer & Gruenwald (1992). 

SV98 provide accurate predictions, as a function of the cluster age, for the total number of 
WR stars and their subtype distribution, the broad stellar emission lines and the luminosities 
and equivalent widths of the two "WR bumps" at A 4650 and A 5808 A . We have used 
these models to derive the age of the ionising population which contains WR stars and whose 
metallicity has been previously derived. 

Regarding the integrated SED from the ionising cluster, the recent models from Leitherer 
et al. (1999; STARBURST99) provide an almost selfconsistent frame to be used in combi- 
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nation with the WR models described above. In fact, they use the same stellar evolution 
models with enhanced mass loss rate, the same atmosphere models to describe the stars in 
the WR phase and cover the range of metallicities and IMF used by SV98. The only appre- 
ciable difference between both sets of models concerns the atmospheres of the main sequence 
stars with initial masses greater than 20 M which in STARBURST99 are represented by 
the plane-paralel Kurucz's models implemented by Lejeune, Cuisiner & Buser (1997). We 
have therefore used the STARBUST99 models in order to fit the emission line spectra for 
our analysed regions. 

3.1 Model fits to the emission line spectra 

In order to predict the emission line spectrum for each analysed region, SEDs predicted by 
the STARBURST99 models for the stellar metallicity closest to the mean gas abundance were 
used as input to the latest version of CLOUDY (Ferland 1999). The models also required the 
gas electron density (assumed to be constant through the nebula) and elemental abundances 
and the ionisation parameter, all of which had been previously derived (see Table 1 and Paper 
I for details). In all cases, we have used as input to the photoionisation models the derived 
values of the ionisation parameter, electron density and epsilon and an inner radius adequate 
to reproduce the derived ionising photons. In all the cases the thickness of the ionised gas 
shell is less than 10% of the total dimensions of the region. This results in a plane-parallel 
geometry. 

We have run models with values of the ionisation parameter inside ± 0.5 dex around the 
derived value and cluster ages with corresponding equivalent temperatures inside ± 2000 K 
around the ionisation temperature estimated from single star photoionisation models. 

The equivalent temperature for a given star cluster has been taken as the stellar ef- 
fective temperature, from Mihalas (1972) model atmospheres, providing the same helium- 
to-hydrogen ionising photon ratio. In all cases the nebulae are assumed to be ionisation 
bounded. Given the run of equivalent temperature with age for ionising clusters,which shows 
an important increase when the most massive stars enter the WR phase, usually more than 
one solution is possible. This can be seen in Figure 1 where we show, for different metal- 
licities 0.2, 0.4 times solar and solar, the cluster equivalent temperature as a function of 
age. 

The best fitting models are then obtained by using an optimisation method that includes 
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the more intense emission line ratios relative to H/?. Our models do not explicity treat the 
physical processes associated with the presence of dust grains. However, since gas-phase 
measured abundances are used as input to the models, the effect of depletion onto dust 
grains is taken into account. The presence of dust grains will affect the thermal equilibrium 
of the nebula, by increasing the electron temperature and producing an enhancement of the 
optical forbidden lines. Nevertheless, according to Shields & Kennicutt (1995), the depletion 
of gas phase coolants is likely to be the most important aspect of dust in terms of its influence 
on the observed spectrum. 

Tables 4 to 14 give, for each region, the predicted line intensities and abundances for 
these models together with the observed ones and their corresponding errors. As can be seen 
most regions show emission line spectra which can be produced by stellar clusters of ages 
between 2 and 3 Myr, although in some cases other solutions are also possible. 

3.2 WR populations 

3.2.1 Region HIS in NGC 628 

The observed WR blue and red 'bump' luminosities relative to H/5 are 0.08 and 0.03 respec- 
tively. According to the models by SV98, maximum values are found for an age of 4Myr 
(0.06 and 0.03 for the blue and red 'bumps' respectively). Moreover, at the age of 4Myr, the 
blue and red 'bumps' equivalent widths are at their maxima, with values around 8.5 A and 
8 A respectively. These values are in excellent agreement with the observed ones: 8.9 A for 
the blue bump and 6.1 A for the red one. Regarding the individual lines, the observed val- 
ues of L(HeII)/L(H/3) and EW(Hell) are 0.04 and 4.8 A respectively. Again, the agreement 
can be considered satisfactory for a single stellar population of about 4Myr (0.02 and 3 A 
respectively for an age of 4.1 Myr). Finally, models predict H/3 equivalent widths of 173 A 
at 4Myr and 159 A at 4.1 Myr, to be compared with the observed value of 140 A. Taking 
into account the global and individual properties of both 'bumps', the agreement between 
the models and the derived quantities is excellent. 

Therefore we can assign an age of 4-4.1 Myr to the ionising stellar population which 
would be responsible for the WR emission features. 
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3.2.2 HII regions in NGC 1232 

From the four HII regions of NGC 1232 analysed in Paper I, three of them present weak blue 
WR 'bumps' (CDT1, CDT3 and CDT4). Special attention must be paid to region CDT1 
showing a metallicity close to solar. For this region the observed WR blue 'bump' luminosity 
relative to H/3 is 0.04, extremely low for a high metallicity HII region. A very low value is also 
found for the blue 'bump' equivalent width, near 2 A. Regarding the stellar Hell luminosity 
and equivalent width, values of 0.015 and 0.8 A are found. According to models by SV98, 
these low values could be produced by an evolved ionising population with an age around 
7 Myr. This population would provide a blue "bump" equivalent width of 2.4 A, consistent 
with the observed value, but would overestimate its luminosity relative to B.f3 by more than 
an order of magnitude (0.55 vs 0.04), while providing an H/? equivalent width of only 4 A 
versus the observed one of 48 A. 

Low values of blue 'bump' luminosities and equivalent widths are also found for early 
ages, at the beginning of the WR phase. A cluster of 2.3 Myr of age provides values of the 
blue 'bump' luminosity 0.02 with respect to H/3 and 5.17 A for its equivalent width, roughly 
consistent with the observations. This population, however, would produce an equivalent 
width of H/3 of 312 A, much higher than observed. 

Similar young ages are found when this sort of analysis is applied to regions CDT3 and 
CDT4. Region CDT3 shows a WR 'bump' luminosity relative to H/? of 0.03. Taking into 
account its derived metallicity (12 + log O/H = 8.56 or Z = 0.008), SV98 predicts this value 
for a cluster with an equivalent width of H/5 of ~ 200A, which is in excellent agreement with 
our measured value of 189 A. This corresponds to an instantaneous burst with age between 
3.1 and 3.5 Myr. Furthermore, the predicted blue 'bump' equivalent width at this age ranges 
from 5.5 to 10 A in good agreement with the observed one (5.9 A). Consistency is also found 
for the weak red 'bump', whose observed values for both the luminosity and equivalent width 
are 0.007 and 2.0 A respectively that should be compared with the predicted ones (0.006- 
0.015 and 2.2-4.2 A respectively). Regarding the Hell line luminosity and equivalent width, 
our observed values are 0.015 and 2.7 A in excellent agreement, again, with the predicted 
ones (0.011-0.007 and 2.15-0.75 A respectively). 

In the case of region CDT4, the WR blue 'bump' luminosity relative to H/5 is 0.08. For 
its derived metallicity (12 + logO/H = 8.37 or Z = 0.006). A value of 0.11 is predicted for an 
instantaneous burst of 4.1 Myr at Z=0.008, which in turn predicts an H/3 equivalent width 
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of 105 A. This value is consistent with the observed one of 138 A. The observed WR blue 
'bump' equivalent width is 10 A in good agreement with the predicted value of 10.8 A at 
this age. The Hell line intensity and equivalent width are again fully consistent with the 
corresponding predicted values (0.03 versus 0.04 and 3.9 A versus 3.7 A, respectively). 

Therefore relatively young single ionising populations between 3 and 4 Myr are required 
to reproduce the observed WR features in regions CDT3 and CDT4. For region CDT1, 
no consistent solution is found. An old cluster of 7 Myr would roughly fit the equivalent 
width of the WR blue 'bump' but would overpredict its luminosity with respect to Hf3 by 
more than an order of magnitude and underpredict the H/3 equivalent width by about the 
same amount. On the other hand, young clusters of about 2.5 Myr predict WR blue 'bump' 
luminosities and equivalent widths in agreement with observations but H/3 equivalent widths 
larger than observed by almost a factor of ten. 

4 DISCUSSION 

The analysis of emission line spectra of GEHR containing WR stars offer an unique oportu- 
nity to test our present ideas about the evolution of star forming regions and check on the 
accuracy of current stellar evolution and atmosphere models. Most work on stellar population 
synthesis of ionised regions assume that the ionising population responsible for the observed 
gas emission belongs to a single star cluster created in a unique star formation episode. 
Therefore the WR emission features can be used to date these star clusters which, in turn, 
should reproduce the observed emission line intensities if the other parameters controling 
the emission line spectra, namely elemental abundances, particle density and ionisation pa- 
rameter, are known. Four of our observed HII regions meet all these requirements. 

4.1 HII regions with observed WR features 

4.1.1 HIS in NGC 628 

H13 in NGC 628, has a mean oxygen abundance of 12+log(0/H)= 8.24, about 0.2 times solar 
if 12+log(0/H)= 8.92 is assumed as the solar value. The models of SV98 for this metallicity 
reproduce the observed WR features for an ionising cluster of 4.1 Myr which produces an 
H/3 equivalent width of 159 A. In this cluster the proportion of WR to O stars is WR/O 
= 5.1xl0~ 2 . Since both SV98 and STARBURST99 models make use of the same stellar 
evolution prescriptions ( see Section 3) we have assumed that a STARBURST99 model of 
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a cluster of metallicity 0.2 times solar which contains the same WR/O star ratio provides 
the spectral energy distribution of the ionising radiation. This model however provides an 
[OIII]A5007 A line emission which is higher than observed by a factor of about 4. Also the 
[SII] and [Nil] lines are down by a factor of about 2 and the equivalent width of H/? is 
somewhat larger than observed. 

On the other hand, our best photoionisation model corresponds to an age of 4.7 Myr 
producing an H/5 equivalent width of 108 A. The most appreciably discrepancy between 
predicted and observed emission line ratios is found for the [OII]A3727 A line which is 
overestimated in the model by a factor of 30% . This line is the one more affected by 
reddening since the rest of the emission lines are measured relative to a nearby hydrogen 
recombination line. The rest of the spectrum is reproduced remarkably well (see Table 4), 
including the equivalent width of H/9 (inside 30 %) and the continuum luminosity at 9000 A 
(inside a factor of two). Clusters younger than 4.7 Myr have a spectral energy distribution 
which results too hard to explain the observations. According to SV98 models a cluster of this 
age would produce WR luminosities and equivalent widths which are lower than observed 
by a factor of two, corresponding to a WR/O ratio of 4.2 xl0~ 2 . Given the uncertainties 
involved in the different models employed, we can conclude that a single star cluster with 
age between 4.0 and 4.7 Myr fit all the observations adequately. The ionising spectrum of 
the cluster of 4.7 Myr is shown in Fig. 2. 

4.1.2 CDT1 in NGC 1232 

Region CDT1 in NGC 1232 shows the highest oxygen abundance of the regions with WR 
features ( 12+log(0/H)= 8.95, about solar if 12+log(O/H) = 8.92). As explained above, 
no consistent solution for a single ionising population is found to reproduce the observed 
WR features. For this region, two star clusters, with 2.3 and 5.7 Myr respectively, provide 
Q(He)/Q(H) ratios nearly identical and corresponding to the stellar effective temperature of 
about 35000 K, in the Mihalas scale, which reproduces the general features of the emission 
line spectrum (see Paper I). Their ionising spectra are shown in Fig. 3. 

The oldest of the two clusters yields emission line ratios in better agreement with the 
observations and an equivalent width of H/3 closer to the observed one (see Table 8). Clusters 
older than this cannot provide the necessary photons to ionise the region. SV98 models for 
this age (5.7 Myr) give however WR features with luminosities relative to H/5 which are a 
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factor of about 10 larger than observed. On the other hand, the young 2.3 Myr cluster is 
more successful in reproducing both the WR blue 'bump' luminosity and equivalent width 
although overpredicts the equivalent width of B./3. These two facts taken together point to 
the presence of a composite population with, at least two clusters: a young one of about 2.5 
Myr containing the WR stars and an older one providing the bulk of the continuum light 
at A 4800 A. We have therefore run a model using as ionising source the combination of the 
spectral energy distributions of two ionising clusters of 2.4 and 7.1 Myr of age calculated with 
the STARBURST99 code in which the youngest of the two provides 10 times the number 
of ionising photons emitted by the oldest. This cluster contains the same WR/O star ratio 
as the SV98 model reproducing the observed WR features. The computed photoionisation 
model is able to reproduce, with a reasonable degree of accuracy, most of the observational 
constraints: the emission line spectrum, the WR features both in luminosity and equivalent 
width, the equivalent width of B./3 and the continuum luminosity at 9000 A. The character- 
istics of the two clusters, together with the corresponding predictions of the photo-ionisation 
model are given in Table 15. The particle density and the elemental abundaces have been 
kept at their derived values given in Table 1. 

4.1.3 CDT3 in NGC 1232 

For region CDT3 in NGC 1232, the SV98 model which best reproduces the observed WR 
features has an age of 3.1 Myr. The ratio of WR to O stars is WR/O=5.4xl0" 2 . A STAR- 
BURST99 model cluster with the same WR/O ratio has 3.3 Myr of age and provides an H/3 
equivalent with of 159 A , close to the observed one (189 A). However, when used as input 
to the photoionisation model, this cluster provides [OIII]A5007 A line emission higher than 
observed by a factor of about 5. This is not unexpected since the equivalent temperature 
corresponding to this cluster is around 38200 K, while that estimated from photoionisation 
models is only 34900 K (see Paper I). Clusters slightly older, also reproducing the observed 
WR features within the errors, have larger WR/O number star ratios (up to 8.3 xl0~ 2 ) and 
therefore their spectral energy distributions result even harder. 

From the analysis of the emission line spectrum, two ionising clusters are found with the 
same Q(He) /Q(H) number photon ratio corresponding to the derived equivalent temperature 
of 34900 K (see Table 10 and Fig. 4). The younger of these clusters however do not contain 



© 2000 RAS, MNRAS 000, 0-?? 



A comprehensive study of reported high metallicity giant HII regions. II. Ionising stellar populations. 

WR stars (WR/O=0.0) and the oldest has WR/O=2.65xl0" 2 , a factor of about 2 lower 
than required to reproduce the WR features. 

From the point of view of ionisation, the combination of two clusters with ages 2.8 and 
4.8 Myr contributing approximately 90 % and 10 % to the total number of ionising photons 
would reproduce the emission line spectrum rather well, except for the [OII]3727 A line which 
results about 40% higher in the model. Also the equivalent width of H/? and the continuum 
light at A9000 A are well reproduced. However, for this population, the luminosity and 
equivalent width of the WR blue 'bump' are lower than observed by factors of about 15 and 
25 respectively. 

4.1.4 CDT4 in NGC 1232 

A similar situation is found for region CDT4. The SV98 model which best reproduces the 
observed WR features has an age of 4.1 Myr. The ratio of WR to O stars is WR/O=9.2xl0~ 2 , 
the highest of all the observed regions. The STARBURST99 model cluster of the same 
age also has the same WR/O ratio and provides an H/? equivalent width of 131 A, that 
should be compared with the observed one of 138 A. When introduced as input in the 
photoionisation code, this cluster reproduces the emission line spectrum satisfactorily, except 
for the [OIII]A5007 A line which is higher than observed by a factor of about 3. Again this 
is related to the hardness of the ionising continuum which provides a Q(He)/Q(H) ratio of 
0.18 which corresponds to an equivalent temperature of 37800 K while the corresponding 
one derived from the emission line spectrum is only 35000 K (see Paper I). 

Again, the analysis of the emission line spectrum indicates two possible star clusters as 
ionising sources whose ages coincide with those found for region CDT3 (see Table 11 and 
Fig. 4). Therefore, regarding the ionisation, the same discussion above applies to this region: 
it is possible to find a combination of ionising clusters which reproduces well the emission 
line spectrum, but predicts WR feature luminosities and equivalent widths well below the 
observed values. 

4.2 HII regions with no WR features observed 

In the rest of the regions no WR features are observed and therefore the estimation of the 
age of the ionising populations must be made with the only help of the analysis of the 
emission line spectra. Out of the two solutions found for region H3 in NGC628, the cluster 
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with 4.5 Myr provides the best fit to the emission line spectrum (Table 5). Most line ratios 
are predicted to better than 10 %. The predicted H/? equivalent width however, is 60 % 
lower than observed. On the other hand, this cluster has a WR/O number star ratio of 
5.8xl0~ 2 , which would provide a WR blue 'bump' luminosity relative to H/? of 0.036 and an 
equivalent width of 3.78 A, that should be observable. The younger cluster, of 2.9 Myr, does 
not have WR stars and, although the fit to the emission line spectrum is slightly worse, it 
does predict adequately both the equivalent width of H/3 and the continuum luminosity at 
A 9000 A. Therefore, a single young ionising cluster of less than 3 Myr, in an evolutionary 
state prior to the appearance of WR stars produces a satisfactory fit to all the observational 
constraints. The spectral energy distribution of the 2.9 Myr cluster is shown in Fig. 2. 

The ionisation of regions H4 and H5 also requires young clusters without WR stars. Clus- 
ters of 2.5 Myr and 2.9 Myr respectively reproduce satisfactorily the emission line spectra, 
the H/3 equivalent widths and the luminosity at 9000 A (see Tables 6 and 7). The cluster in 
H5 seems to be slightly older than that in H4 and the ionisation parameter in this region is 
lower by a factor of about 2. The SEDs of the two ionising clusters is shown in Fig. 5. 

The region in NGC 1232 not showing WR features, CDT2, also allows for two possible 
ionising clusters (see Table 9 and Fig. 6). According to the STARBURST99 models, the 
oldest one with 4.5 Myr has a WR/O star ratio of 0.1 which, in SV98 models, gives WR 
blue 'bump' luminosity and equivalent width of 0.1 and 6.8 A respectively. On the other 
hand, the youngest cluster with 3 Myr, has a WR/O star ratio of only 0.01 and provides WR 
blue 'bump' luminosity relative to H/3 and equivalent width of 0.01 and 2.34 A respectively. 
Based on the fact that no WR features are observed, we consider this latter cluster as a 
better candidate for the ionisation of the region. Nevertheless, it does not seem to be the 
only population present since the low observed value of the H/5 equivalent width points to 
an underlying population that would contribute 70 % to the luminosity at H/3. The WR blue 
'bump' equivalent width would be consequently lowered to only 0.7 A which is consistent 
with the lack of detection of the feature. 

The other region of the sample showing a high metallicity (about twice solar) is CDT1 in 
NGC 1637. Two ionising clusters with Z=0.04 are able to reproduce the observed emission 
line spectrum to the same degree of accuracy and they produce H/? equivalent widths higher 
and lower than observed respectively (see Table 12 and Fig. 7)). Also the observed A9000 
A luminosity is intermediate to that predicted by the two clusters. Therefore, the best 
solution would seem to consist of a combination of the two. We should however keep in 
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mind that, according to models, both clusters contain WR stars with WR/O ratios between 
0.14 and 14. The oldest cluster would provide large values of both blue 'bump' luminosity 
and equivalent width (0.7 and 13 A) while the younger one would provide about the same 
equivalent width but a luminosity lower by a factor of about 5 (0.14). Therefore, based on 
the fact that no WR features are observed we feel that the most plausible solution would 
be the combination of a young ionising cluster providing the bulk of the ionisation and an 
underlying population which would lower the equivalent width of H./3 to the observed values. 
This underlying population should contribute 50% of the continuum luminosity at A 4861 A. 

Finally, not much can be said about the two regions in NGC 925. Both of them show 
a metallicity of about 12+log(0/H)~8.5, although, due to the non detection of suitable 
temperature sensitive lines, the uncertainty in the derived oxygen abundance is high (± 0.20). 
In the case of CDT1, the emission line spectrum can be adequately reproduced by relatively 
old ionising clusters of slightly different age and metallicity (5.8 Myr and 12+log(0/H) = 
8.52 or 5.9 Myr and 12+log(0/H) = 8.62) with the best solution probably being in between 
(see Table 13). There is however some room for contribution by a non ionising cluster since 
the luminosity at A9000 A is predicted to be lower than observed by a factor of about 3. 

In the case of CDT4, the 5.9 Myr cluster predicts an equivalent width of B./3 lower than 
observed by a factor of about two, while a young cluster without WR stars (2.8 Myr) gives 
a value much larger than observed (Table 14). It is however possible to find a satisfactory 
solution by considering the combination of two clusters of 2.8 and 5.9 Myr with this former 
one providing 60% of the total number of ionising photons. This solution is also listed in 
Table 14. 

4.3 Global analysis 

In summary, out of the 11 analysed regions, only the four in NGC 628 seem to be consistent 
with the presence of single ionising star clusters, all of them show relatively low abundances. 
For the rest of the regions, at least two populations are needed: a young one, with age between 
2.5 and 3 Myr, which provides most of (or all) the ionising photons, and an older one which 
contributes a substantial fraction of the continuum luminosity. This second population also 
contains an important fraction of the mass of the composite cluster. The difference in age 
between the two populations is not very large, 3-10 Myr. 

Composite populations for HII regions have been found in previous work by Mayya & 
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Prabhu (1996) and by Diaz et al. (2000b) for disc and circumnuclear objects respectively 
from broad band and Ha photometry A spectrophotometric analysis of circumnuclear HII 
regions in NGC 3310 (Pastoriza et al. 1993) and NGC 7714 (Garcia- Vargas et al. 1997) also 
reveals the presence of composite populations and even in the young HII region 30 Dor in 
the LMC a spread in age for the stellar populations has been found (Selman et al. 1999). 

This scenario, is able to reproduce both the emission line spectrum and the WR features 
of the high metallicity region CDT1 in NGC 1232. It is nevertheless not sufficient to explain 
the observations in the two intermediate metallicity HII regions in NGC 1232 with WR 
features: CDT3 and CDT4. For these two regions, no consistent solution can be found. While 
the models by SV98 consistently provide excellent fittings to the observed WR properties, 
the evolutionary models by STARBURST99 result too hard to explain the observed emission 
line spectra. 

Several solutions can be suggested such as the inadequacy of the stellar evolutionary 
tracks, uncertainties in the computation of stellar atmosphere models or the everlasting 
question of temperature fluctuations within these regions that would underestimate the true 
abundances. 

The fact that the WR features are adequately reproduced by SV98 models seem to imply 
that the evolutionary tracks are able to predict the right relative numbers of WR/O stars, 
and their different subtypes at the derived abundances. These relative numbers, combined 
with the observed emission line luminosities of the individual WR stars, and the predicted 
continuum energy distribution of the ionising population predict emission line intensities 
and equivalent widths of the WR stars that are in excellent agreement with observations. 
This can be seen in Fig 9. Upper and lower panels respectively show the predicted emission 
line intesities and equivalent widths of the WR blue 'bump' for different metallicities: 0.2 
solar, 0.4 solar and solar together with the data corresponding to the four observed objects. 
Solid symbols correspond to the observed data. In the case of CDT1 in NGC 1232, for which 
a composite population is found, two open symbols are shown, one corresponding to the 
predictions for a single ionising cluster, and another one including the contribution to the 
H/3 and continuum luminosity of the cluster not containing WR stars. For regions CDT3 
and CDT4 observations and predictions are nearly identical for the respective clusters and 
therefore only one symbol is shown. As we can see, the agreement for the three regions in 
NGC 1232 is surprisingly good. For region H13 in NGC 628 predicted and observed values 
(solid and open symbols) agree inside a factor of about two. We therefore conclude that, 
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regarding the WR star populations, the currently assumed stellar evolution, including a 
high mass loss, is adequate to match the observations, at least for intermediate and high 
metallicities. The worse agreement found for region H13 in NGC 628, whose WR feature 
luminosity and equivalent width seem to indicate a slightly higher metallicity, might in 
principle be attributed to the presence of temperature fluctuations, which do not need to be 
invoked in the other cases. 

The third possible source of disagreement concerns the existing uncertainties in modelling 
the atmospheres of high mass stars of O and WR types. Our results seem to point to an 
overestimate of the hardness of the spectral energy distribution of the stars in the WR phase 
which is more apparent for intermediate to high metallicities. The same sort of effect has been 
found by DCTGOO for region 74C in NGC 4258 and Esteban et al. (1993) for the galactic 
WR nebula Ml-67. These latter authors, from stellar and nebular spectroscopic analyses, 
concluded that lower temperatures were required from the photoionisation models for late 
type WN (WNL) stars. In their study they used the unblanketed WR models of Schmutz et 
al. (1992). A subsequent reanalysis of this region has been made by Crowther et al. (1999) 
using blanketed model atmospheres. In this case the resulting ionising spectrum is much 
softer and reproduces better the observations although some important discrepancies, like a 
negligible Hel A 5876 A line intensity, which is nonetheless observed, and the overprediction 
of [Oil] A 3727 A and [Nil] A 6584 A by factors of 3 and 2 respectively, still remain. 

It should however be recalled that all results discussed above are based on ionisation 
bounded models. An interpretation of the emission line spectra of HII region with WR 
features on the basis of matter bounded models has been presented in Castellanos, Diaz & 
Tenorio-Tagle (2002). Models of this kind have been found to provide excellent fittings to 
the observations of H13 in NGC 628, CDT3 in NGC 1232 and 74C in NGC 4258 using both 
SV98 and STARBURST99 models. These models would point to an important leakage of 
ionising photons depending on both the metallicity and evolutionary state of the region. 

Certainly, the inclusion of line blanketing in the model atmospheres of WR stars, which 
is already included in the CoStar (Schaerer & de Koter 1997) and Lejeune et al. (1997) 
models is imperative to interpret correctly the emission line spectra observed in HII regions. 
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5 SUMMARY AND CONCLUSIONS 

We have analysed the stellar populations of eleven Hn regions in the galaxies NGC 628, 
NGC 1232, NGC 925 and NGC 1637 using spectrophotometric observations between 3500 
and 9700 A. We have derived the physical properties of the regions and their corresponding 
ionising clusters: filling factor, mass of ionised gas and mass of ionising stars. Most of the 
regions have small ionising clusters with masses in the range 1500 to 30000 solar masses. 
The exceptions are the three supergiant HII regions in NGC 1232, CDT1, CDT3 and CDT4 
with masses greater than 100,000 solar masses. These values constitute in fact lower limits 
since the regions are assumed to be ionisation bounded and the presence of dust has not 
been taken into account. 

WR features have been observed in the four supergiant HII regions, H13 in NGC 628 
and CDT1, CDT3, CDT4 in NGC 1232. A very detailed modelling has been carried out for 
these regions by using two different sets of models: those of SV98 for WR star populations 
and those of Leitherer et al. (1999) for ionising populations, to try to reproduce simulta- 
neously both the WR features and the emission line spectrum of each region. In all cases, 
the agreement between both predicted and observed values for the WR luminosities and 
equivalent widths is excellent. This fact seems to indicate that the stellar evolution assumed 
by the models is predicting the right WR/O ratios at the different metallicities involved. 

The ages of the populations containing WR stars are found to be between 2.3 and 4.1 
Myr. This population is sufficient to explain the emission line spectrum of the low metallicity 
region H13 in NGC 628, taking into account the uncertainties involved in both observations 
and model computations. 

This is not the case for the rest of the regions. For CDT1 in NGC 1232, a second ionising 
population is required to simultaneously reproduce both the WR features and the emission 
line spectrum. This second population is older than the previous one by about 5 Myr and 
contributes about 10% of the ionising photons and 5 times the continuum luminosity at 
9000 A. Composite populations are also found for half of the regions without WR features, 
although in this case, the result is based only on the emission line spectrum analysis. 

For the other two regions containing WR stars, CDT3 and CDT4 in NGC 1232, no 
consistent solution is found, since the population containing WR stars produces a spectral 
energy distribution which is too hard to explain the emission of the gas. Several solutions are 
proposed to solve this problem. Both a reduction in the number of high energy photons of 
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the ionising clusters, which could arise naturally with the inclusion of blanketing by metallic 
lines in the atmospheres of the WR stars, and/or the assumption that the HII regions be 
matter bounded can effectively solve the problem. The latter hypothesis has been tested 
successfully (Castellanos et al. 2002) while the former one still remains to be explored. 
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Table 1 . Summary of derived gas physical conditions of the analysed HII regions 
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Table 2. Physical properties of the observed HII regions 
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1.77 


13.0 






10586 


1.5 



Table 3. WR feature intensities and equivalent widths in the observed HII regions. 

Region L(WR)/H/3 EW(WR)(A ) L(HeII)/H/3 EW(HeII)(A ) 

Region H13 (NGC 628) 0.08 8.9 0.04 4.8 

Region CDT1 (NGC 1232) 0.04 2.0 0.015 0.8 

Region CDT3 (NGC 1232) 0.03 5.9 0.015 2.7 

Region CDT4 (NGC 1232) 0.08 10 0.03 3.9 
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Table 4. Evolutionary models for region H13 in NGC 628 
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Table 5. Evolutionary models for region H3 in NGC 628 
Parameter Observations 2.9 Myr 



Q(He)/Q(H) 




0.13 


Q(He+)/Q(He) 




0.00 


WR/O 




0.0 


<log U> 


-2.92 ± 0.15 


-2.89 


n e 


<40 


10 


12 + log(0/H) 


8.23 ± 0.15 


8.23 


log(S/0) 


-1.85 ± 0.11 


-1.73 


log(N/0) 


-1.14 ± 0.10 


-1.11 


3727 [Oil] 
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5007 [OIII] 
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11 
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Table 6. Evolutionary models for region H4 in NGC 628 
Parameter Observations (H4) 2.5 Myr 



Q(He)/Q(H) 0.18 

Q(He+)/Q(He) 0.00 

WR/O 0.0 

<log U> -2.95 ± 0.15 -2.95 

lie <40 10 

12 + log(0/H) 8.31 ± 0.15 8.18 

log(S/0) -1.80 ± 0.11 -1.71 

log(N/0) -0.97 ± 0.10 -1.06 

3727 [Oil] 4432 

5007 [OIII] 2177 ±15 2000 

6584 [Nil] 527 ± 20 500 

6717 [SII] 308 ± 20 245 

9069 [SIII] 147 ± 12 222 

6312 [SIII] 10 ± 2 16 

t(0 2 +) 0.94 1.09 

t(S 2 +) 0.95 ± 0.05 1.11 

EW(H/3)(A) 220 290 

L9000 (10 34 erg s" 1 !" 1 ) 2.9 4.2 



Table 7. Evolutionary models for region H5 in NGC 628 
Parameter Observations (H5) 2.95 Myr 



Q(He)/Q(H) 
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0.79 


0.97 
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0.82 ± 0.06 


1.0 
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222 
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3.4 
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Table 8. Evolutionary models for region CDT1 in NGC 1232 
Parameter Observations 2.3 Myr 5.7 Myr 
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Table 9. Evolutionary models for region 


CDT2 in NGC 1232 


Parameter 


Observations 


3.0 Myr 


4.5 Myr 


Q(He)/Q(H) 




0.155 


0.130 


Q(He+)/Q(He) 




5.4xl0~ 4 


0.062 


WR/O 




0.01 


0.1 


<log U> 


-2.95 ± 0.15 


-3.05 


-3.10 


n e 


<40 


10 


10 


12 + log(0/H) 


8.61 ± 0.15 


8.60 


8.55 


log(S/0) 


-1.81 


-1.71 


-1.71 


log(N/0) 


-1.15 


-1.15 


-1.15 


3727 [Oil] 


4180 ± 220 


4313 


4026 


5007 [OKI] 


1350 ± 20 


1375 


1334 


6584 [Nil] 


717 ± 30 


710 


694 


6717 [SII] 


537 ± 15 


517 
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191 ± 15 
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0.81 
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Table 10. E volution ary 


models for region 


CDT3 in 


NGC 1232 


Parameter 


Observations 


2.8 Myi 


4.8 Myr 


Q(Hc)/Q(H) 




0.088 


0.082 


Q(He+)/Q(He) 




5.1x10" 


6 0.25 


WR/O 




0.0 


0.027 


<log U> 


-2.72 ± 0.10 


-2.80 


-2.82 


n e 


223 


230 


230 


12 + log(0/H) 


8.56 ± 0.14 


8.58 


8.55 


log(S/0) 


-1.70 ± 0.07 


-1.70 


-1.71 


log(N/0) 


-0.96 ± 0.04 


-1.00 


-0.96 


3727 [Oil] 


3180 ± 150 


4471 


2868 


5007 [OKI] 


842 ± 13 


823 


1972 


6584 [Nil] 


892 ± 30 


894 


890 


5755 [Nil] 


8 ± 2 


10 


9 


6717 [SII] 


333 ± 10 


301 


329 
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48 ± 6 


10 


43 
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229 ± 10 


340 


271 


6312 [SIII] 


8 ± 2 


13 


10 


t(0 2 +) 


0.69 


0.82 


0.81 


t(S 2 +) 


0.74 ± 0.05 


0.85 
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0.89 
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<log U> 
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n e 
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100 
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8.37 ± 0.12 


8.45 
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0.90 


0.87 
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0.90 ± 0.06 


0.94 
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Table 12. Evolutionary 


models for region 


CDT1 in 


NOC 1637 


Parameter 


Observations 


2.8 My) 


; 5.9 Myi 


Q(He)/Q(H) 




0.082 


0.064 


Q(He+)/Q(He) 




1.26x10" 


" 4 1.20x10' 


WR/O 




0.24 


>0.9 


<log U> 


-3.30 ± 0.20 


-3.25 


-3.15 


n e 


100 


100 


100 


12 + log(0/H) 


9.10 


9.20 


9.10 


log(S/0) 


-1.97 


-1.90 


-1.80 


log(N/0) 


-0.88 


-0.87 


-0.87 


3727 [Oil] 


1170 ± 60 


1116 


1031 


5007 [OKI] 


84 ± 10 


90 


90 


6584 [Nil] 


1037 ± 60 


1212 


1039 


6717 [SII] 


335 ± 20 


359 


379 


9069 [sin] 


99 ± 9 


131 


141 




0.40 


0.46 


0.46 


EW(H/3)(A) 


74 
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22 


L9000 (10 35 erg s^A" 1 ) 


2.6 


1.1 
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CDT1 in 


NOC 925 


Parameter 


Observations 


5.8 Myi 


: 5.9 My] 


Q(He)/Q(H) 




0.124 


0.076 


Q(He+)/Q(He) 




7.58x10" 


- 4 2.86x10" 


WR/O 




0.023 


0.022 


<log U> 


-3.00 ± 0.30 


-3.15 
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n e 
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log(N/0) 
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3727 [Oil] 
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849 ± 14 
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6584 [Nil] 
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Table 14. Evolutionary models for region CDT4 in NGC 925 
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Table 15. Ionising populations for region CDT1 in NGC 1232 



Parameter 


2.4 Myr 


7.1 Myr 


Composite 


Q(H) (10 50 s" 1 ) 


13.8 


1.38 


15.2 


Q(He)(10 48 s- 1 ) 


140 


0.17 


140 


L(H/3)(10 37 erg s" 1 


66.1 


6.61 


72.7 


M*(10 4 MO) 


2.73 


15.3 


18.0 
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0.0 


0.026 


EW(WR)(A) 


7.3 
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Figure 1. The equivalent temperature, defined as that of the Mihalas model atmosphere which provides the same ratio of 
helium to hydrogen ionising photons, versus age for single ionising clusters of different mctallicity. 
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Figure 2. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectra of regions H13 and 
H3 in NGC 628 
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Figure 3. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectrum of region CDT1 
in NGC 1232 
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Figure 4. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectra of regions CDT3 
and CDT4 in NGC 1232 
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Figure 5. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectra of regions H4 and 
H5 in NGC 628 
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Figure 6. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectrum of region CDT2 
in NGC 1232 
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Figure 7. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectrum of region CDT1 
in NGC 1637 
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Figure 8. Spectral energy distributions of the ionising clusters able to reproduce the emission line spectra of regions CDT1 
and CDT4 in NGC 925 
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Figure 9. Relative intensity (left) and equivalent width (right) of the WR blue 'bump' versus H/3 equivalent width for SV98 
models of three different metallicities as labelled. The data are shown as solid and open symbols as explained in the text. 
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